Vortex Deconfinement in the XY Model with a Magnetic Field 
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We study vortex unbinding for tlie classical two-dimensional XY model in a magnetic field on square 
and triangular lattices. A renormalization group analysis combined with duality in the model shows 
that at high temperature and high field, the vortices unbind as the magnetic field is lowered in a two- 
step process: strings of overturned spins first proliferate and then vortices unbind. The transitions 
are highly continuous but are not of the Kosterlitz-Thouless type. The unbound vortex fixed point 
is shown to inherit properties of the underlying lattice, in particular containing a set of nodal lines 
that reflect the lattice symmetry. 
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I. INTRODUCTION 

Topological defects play a crucial role in a host of phe- 
nomena in condensed matter systems^. Among these are 
Kosterlitz-Thouless (KT) transitions^, in which vortex- 
like defects in a two dimensional system unbind above 
some critical temperature. Kosterlitz-Thouless transi- 
tions are important for understanding disordering in pla- 
nar (XY) magnets, destruction of superfluidity and su- 
perconductivity in thin films, melting of two dimensional 
crystals, and transitions among different liquid crystal 
states^. Vortex unbinding is related to the roughening 
transition, in which the interface between two materials 
fluctuates without bound above a critical temperature'^. 
They are important for understanding ground state prop- 
erties of 1-1-1 dimensional quantum systems (Luttinger 
liquids)^. They have also recently found relevance in un- 
derstanding the states of stripe systems in high tempera- 
ture superconductors and quantum Hall systems^'^. The 
XY fcrromagnct is often considered the paradigm of all 
these systems, largely because of the (somewhat mislead- 
ing) simplicity of its Hamiltonian. 

Much less attention has been paid to what becomes 
of the vortex unbinding transition in the presence of a 
symmetry-breaking field, such as a magnetic field tending 
to align the XY spins along some particular direction, in 
spite of the fact that most of the systems discussed above 
have realizations in which such symmetry-breaking fields 
are present. One major reason for this is that the stan- 
dard KT analysis breaks down under the most interesting 
circumstances when a symmetry breaking field is present. 
Usually, one assesses whether vortices may be incorpo- 
rated into a simple theoretical description of a system by 
renormalization of the parameters of the Hamiltonian. 
This is possible when the vortices are bound tightly into 
pairs of vanishing net vorticity. In the renormalization 
group (RG) approach to the KT transition, the paired 
state is signaled by a flow in which the vortex fugacity, 
Q-Ea/ksT ^ y^2^ with Ec the core energy of a vortex, 
vanishes. The interpretation of this is that an effective 
long- wavelength theory of the system may be constructed 
with vanishing vortex density (y = 0); the vortex degrees 



of freedom are irrelevant in the renormalization group 
sense^. 

For small values of y and above the Kosterlitz-Thouless 
transition temperature Tkt, one finds that y grows 
rather than shrinks, signaling the relevance of vortices 
in the system, so that a theory without vortices cannot 
be qualitatively correct. This sharp change in behavior 
is accompanied by an essential singularity in the free- 
energy, so that vortex unbinding constitutes a continu- 
ous phase transition. Because the RG flows contain a set 
of simple attractive fixed points (the Hamiltonians with 
y = to which the low temperature microscopic Hamil- 
tonians fiow under the RG), thermodynamic properties 
of the transition may be computed^'^. 

For many systems with symmetry- breaking fields, the 
existence of this set of stable fixed points is lost®. A 
prototype of this - and the system we focus on in this 
paper - is the XY system with a magnetic field. It is 
defined by the Hamiltonian 



nxY/kBT = -K J2 



Or') - h^COsOr. (1) 



<r,r'> 



In Eq. 1, 6r represents the angle of a planar spin at 
location r, K is an exchange coupling between nearest 
neighbor spins, and h is the magnitude of an effective 
magnetic field tending to align the spins along the x 
axis. For convenience from now on we will choose our 
units of energy so that ksT =1. A perturbative RG 
analysis® when applied to this system shows that either 
the symmetry-breaking term in the Hamiltonian or the 
vortex fugacity, or both, are relevant even when they are 
initially very small. The absence of a stable fixed point 
means this perturbative analysis cannot provide us with 
a simple Hamiltonian to describe the state of the system. 
This has been interpreted to mean that there is no phase 
transition, and presumably no possibility of vortex un- 
binding, when a symmetry-breaking field is present (see, 
for example. Chap. 2 of Ref. 1). In this paper, we will 
describe a new analysis that suggests that while the first 
part of this interpretation (no phase transition) is at some 
level true, the second part (no vortex unbinding) is not: 
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we will develop a RG description that allows for vortex 
unbinding, but with a fixed point structure that remark- 
ably avoids the free energy singularity usually associated 
with a phase transition. Some of the results discussed 
here were described elsewhere'^. In this paper we provide 
details of those calculations, and extend them to show 
how the underlying lattice symmetry may be incorpo- 
rated into the model. 

An important perspective on why unbinding in the XY 
ferromagnet is different in the presence of a magnetic field 
than without it comes from examining the interaction of a 
single vortex-antivortex pair. Without symmetry break- 
ing, this is well known to be a slow, logarithmic poten- 
tial with respect to pair separation. This arises because 
the lowest energy spin configuration, subject to the con- 
straint that the system contains a vortex and an antivor- 
tex, approaches the ferromagnetic ground state slowly as 
one moves away from the pair. Such a configuration is 
energetically expensive for h > since it involves many 
spins tilted away from the direction specified by the field; 
the system instead forms a string of overturned spins con- 
necting the vortex-antivortex pair-'^'^. The rotation of the 
spins through 27r as one moves through the string is es- 
sentially identical to a kink soliton of the sine-Gordon 
model^^, carrying a ch aracteristic energy per unit length 
proportional to VKh. Thus, the potential energy for a 
pair of vortices increases linearly with separation. Such 
linear potentials arise in the context of strong interac- 
tions, and leads to confinement^^, a situation in which 
individual particles (quarks for strong interactions, vor- 
tices in the XY model) do not appear in isolation. From 
this analogy one might conclude there is no unbound vor- 
tex state of this system. 

The difficulty with this argument is that it ignores 
the statistical mechanics of the string connecting the 
vortices. At high temperatures, one must include en- 
tropy in the free energy of the string connecting a sin- 
gle vortex-antivortex pair. The number of configurations 
for a string grows exponentially^^ with its length L, so 
that at high temperature a string cannot bind a single 
vortex-antivortex pair. Moreover, the string tension its 
energy per unit length - is renormalized downward by 
small breaks (vortex-antivortex pairs), which may occur 
along its length. Finally, at high temperature the XY 
system will be flooded with vortices and antivortices, so 
that the distance between neighboring pairs will be of the 
same order as the pair sizes themselves. This situation 
is highly reminiscent of what is expected to occur in a 
quark-gluon plasma, in which quarks are deconfined^^. 
Deconfinement at high temperature is clearly also a pos- 
sibility for vortices in the XY system. To decide if and 
when this occurs, we will develop a model that deals di- 
rectly with the strings, and contains a parameter which 
we will identify with the string tension. This string ten- 
sion can be driven to zero when renormalized by fluctua- 
tions - signaling an unbinding transition for the vortices. 

Before beginning our technical discussion, we summa- 
rize our results. We develop an effective model for this 



system, focusing on a representation in which the string 
degrees of freedom are explicit. The model has a dual 
form, which may be interpreted as a solid-on-solid model 
with screw dislocations. The resulting low energy exci- 
tations c;ontain both open and closed domain walls. The 
endpoints of the domain walls - the screw dislocations - 
are degrees of freedom dual to the vortices. We initially 
focus on the states of these dislocations, and then use the 
duality to draw conclusions about the allowed states of 
the vortices. The perturbative RG analysis we follow is 
valid in the large h and small y limit. It generates a pa- 
rameter p, which represents the energy per unit distance 
for separating two endpoints of a single open domain wall, 
keeping the length of the domain wall fixed. The RG 
equations indicate that p may flow to zero, or may flow 
to a flnite value; we interpret the latter as a bound dislo- 
cation state. The p = fixed point is accessible when the 
domain walls are rough, i.e., when they are unbounded in 
size and percolate through the system. Since arbitrarily 
large domain walls may be cut open with endpoints as far 
apart as we like without extra energy cost at the p = 
fixed point, we identify this as an unbound dislocation 
phase. 

In the RG fiows, the p = fixed point appears at the 
end of a line of fixed points with p> 0, so that the flows 
accumulate at the p = for the unbound phase. The 
string tension p grows continuously from zero when the 
microscopic parameters of the Hamiltonian cross into the 
values corresponding to the bound dislocation phase. For 
fixed K and Ec, bound dislocations occur for h smaller 
than some critical value. 

An important observation is that there is no relevant 
direction in the parameter space of the Hamiltonian lead- 
ing away from the unbound dislocation fixed point. This 
is a remarkable result: in the absence of a relevant direc- 
tion, there is no mechanism by which the effective free 
energy can accumulate a singularity in the RG as one 
integrates out short distance scales. With no free energy 
singularity, one does not expect to find singularities in 
any thermodynamic quantities for the system. Neverthe- 
less, we can sharply distinguish the bound and unbound 
vortex phases, for example by measuring the diffusion 
constant for vortices, or (equivalently) measuring fluctu- 
ations in the vortex dipole moment of the system. Recent 
Langevin dynamics simulations^^ focusing on the latter 
quantity have confirmed the basic results of the RG stud- 
ies presented here; we will discuss these in more detail 
below. The important lesson at this point is that vortex 
deconfinement in this system is not a phase transition in 
the usual sense. One can distinguish the bound and un- 
bound vortex phases through transport properties, cor- 
relation functions, or statistics of specific fluctuations in 
the system. But they cannot be distinguished by quali- 
tative differences in quantities that may be expressed as 
derivatives of the free energy; i.e., thermodynamic quan- 
tities. The possibility of such unusual transitions has 
been noted in rigorous studies of phase transitions^^, but 
to our knowledge this is the first concrete example of such 
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behavior. 

The unbound dislocation fixed point, we shall see, has 
a unique signature: it contains one or more lines of zero 
energy modes in the Brillouin zone. Our earlier analysis^ 
focused on the case of a single such nodal line. This 
breaks the discrete rotational symmetry of the lattice, 
and we will see it arises from using a particular choice of 
gauge which leads to approximations breaking the sym- 
metry. We will show that the underlying symmetry of 
the lattice may be retained without changing any of the 
qualitative results. In the case of a square lattice, the 
fixed point contains two nodal lines, and for a triangular 
lattice it contains three. The number of nodal lines is 
determined by the number of independent directions a 
domain wall may exit from a dislocation; it is essentially 
the number of distinguishably different types of dislo- 
cations the system supports. We note that these fixed 
points might be regarded as a classical analog of simi- 
lar Hamiltonians arising in the context of quantum spin 
models with ring exchange 

To understand what all this implies about vortex un- 
binding, we must re-express the partition function in its 
dual form. The dual Hamiltonian we will see is identi- 
cal to the dislocation Hamiltonian, provided we exchange 
Ec ^ l/2h and substitute K l/Air'^K. Since the un- 
bound dislocation phase occurs for large Ec and small 
l/2/i, we expect unbound vortices can be found at large 
l/2h and small Ec- That an unbound vortex phase does 
not occur for large values of Ec shows how different the 
transition is from the KT phenomenology, and also is con- 
sistent with the failure of perturbative RG calculations^ 
with small y = e~^° (i.e., large Ec) to capture the transi- 
tion. Since dislocation unbinding occurs at a finite value 
of l/2/i, dislocation unbinding and vortex unbinding can- 
not be the same transition. This tells us there must be an 
intermediate phase in which both dislocations and vor- 
tices are bound. The nature of this phase follows from 
the observation that domain walls arc rough even in the 
bound dislocation phase. At long wavelengths, we can 
ignore the dislocations, so one is effectively in the rough 
(high temperature) phase of a solid-on-solid model. This 
is dual to a bound (low temperature) phase of logarith- 
mically interacting vortices^, so we identify the large Ec, 
l/2h limit of the parameters as a logarithmically bound 
dislocation phase. From the duality, we expect this im- 
plies the vortices will also be logarithmically bound in 
this region of parameters. Physically, we can understand 
this by recognizing the domain walls in the dislocation 
representation are dual to the strings in the vortex repre- 
sentation, and both are proliferated for these parameters. 
Proliferated strings do not linearly confine the vortices, 
but they can (and do) upwardly renormalize the loga- 
rithmic interactions between them, enough so that the 
vortices may be bound together even if T > Tkt- 

Finally, one may ask what the unbound dislocation 
phase looks like when expressed in terms of the vortex 
degrees of freedom. Since this occurs for very large values 
of h and Ec, the bare string tension between vortices is 




l/2h 

FIG. 1. Schematic phase diagram for XY model in a mag- 
netic field for K « l/27r. Upper left corner represents linearly 
confined phase (unbound screw dislocations in dual represen- 
tation). Middle phase contains proliferated loops in both 
descriptions and vortices (dislocations) are logarithmically 
bound. Lower right corner contains unbound vortices (lin- 
early confined dislocations). 

very large and the vortices are dilute. In this case it 

is clear that the vortices are linearly confined. Thus we 
identify the highest temperature phase of the dislocations 
- a deconfined phase - with the lowest temperature phase 
of the vortices - a linearly confined phase. 

All these considerations suggest that for large Ec 
and/or small h, there are three phases for the vortices: a 
linearly confined phase, a logarithmically confined phase, 
and a deconfined phase. This is illustrated in Fig. 1, 
which is the simplest phase diagram one may draw con- 
sistent with the perturbative RG analysis. 

We are led to a picture in which vortices deconfine in 
a two-step process: as fluctuations increase, the strings 
connecting them proliferate, but vortex-antivortex pairs 
remain bound by a residual logarithmic attraction; at 
still higher levels of fluctuations the closed string loops 
may break open, and the vortices deconfine. We note 
that simulations strongly support the existence of two 
different bound phases, as well as the unbound vortex 
phase^^. 

This article is organized as follows. We begin Section II 
by expressing the system as a Villain model, and then for- 
mulate the model as a continuum theory. In Section III 
wc develop the RG calculations and show how the un- 
usual fixed point structure emerges. Section IV discusses 
more careful formulations that respect the lattice symme- 
try for square and triangular lattices. Section V contains 
discussion of the numerical integration of the RG equa- 
tions, and we conclude with a summary in Section VI. 
Two Appendices provide further details of the calcula- 
tions. 
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II. DERIVATION AND INTERPRETATION OF 
THE MODEL 

A. Representation in Terms of Domain Walls and 
Screw Dislocations 

In this section, we present the derivation of the model 
we will be considering. Many of the steps we take may be 
found in Ref. 8, and indeed our final result is equivalent 
to theirs, although the representation on which we finally 
perform our RG analysis is considerably different, leading 
to different results. In any case, for completeness we 
present the derivation in its entirety. 

The partition function for the XY model may be ex- 
pressed as ZxY = J V9e~^^^ . Because of the cosines 
appearing in the exponent, it is very difficult to make 
progress working directly with this Hamiltonian. Since 
we are interested in phases and long wavelength proper- 
ties of the system, we can replace the model with any 
other that contains the same symmetries and retain the 
correct physics. A particularly useful model to adopt is 
the Villain model^'^, in which one makes the replacement 



exp — 27rZ^^Tr^r — i (^r — Or')SrT' 

r (r,r') 

-iJ2^rTr], (4) 



where Srr' and are the integer fields arising from the 
Poisson resummations. The continuous fields in Eq. 4 

may now all be integrated out. Wc focus first on the 
angular variables 9, for which the integration produces a 
product of (5-functions. These may be conveniently writ- 
ten if we re-express the bond integer field S^r' as a vector 
field S(r), with components Si{r) = Sr,T+Ai, where Aj 
are vectors denoting nearest neighbor bonds in the "posi- 
tive" direction (For example, on a square lattice one may 
take Ai = aoy, A2 — aox.). It is important to note that 
the dimensionality of S is not set by the (two-) dimen- 
sionality of the lattice, but rather is given by (half) the 
coordination number; for example, in the triangular lat- 
tice discussed below, S is a thrc;e-dimcnsional vector field. 
Whatever its dimensionality, we can define a divergence, 

V-S(r)^^[5,(r)-5,(r-A,)]. 



„Ccos0 



E 



-C{0-2Trmf/2 



(2) The second line of Eq. 4 can now be rearranged to read 



whenever a cosine appears in the exponent, m above is 
an integer variable, and the important observation is that 
the weight as a function of 9 has period 2tt on both sides 
of Eq. 2. The replacement of a cosine in the exponent 
with a quadratic form allows us to make some progress. 
Applying this to Zxy, we are led to the replacement 



2xY ZvM 



K 



~y E (^r - ^r' - 27rmr,r')^ 



'-X9,-2Tmrf 



(3) 



The functional integral over the 9 field may be performed 

with the help of the Poisson resuinmation formula. This 
states that for a function g{m) summed over an integer 
field m, 



E 9{m)= 



m— — 00 



n— — 00 



d(j) g{(j))e 



Applying this to both integer fields in Eq. 3 and shifting 
variables, we arrive at the expression 



2vM = J2Yl 1^9 Vcl> 

S„, Tr. 

X exp [ -2'k'^K - 27r^/i^ Vr 



(r,r') 
2TTi Y S, 



rr' V^riT' 



X exp 



-27rj Y ^rV-r -iY^r[Tr-V ■ S(r)] 



from which it is now clear that the angular integration 
yields, up to unimportant factors of 27r, Hr ^Tr,v s(r)- 
Substituting this into Eq. 4, performing the remaining 
continuous integrals and using the 6 functions to elimi- 
nate the T sum, we arrive at the remarkably simple form 
(ignoring unimportant prefactors) Zvm = J2s e~^^** , 
with 

^VM = ^ E ISWI^ + ^ E(V • S(r))'- (5) 



In the limit /i — > 0, the configurations entering the par- 
tition sum must satisfy V • S(r) = 0, which implics"'^* one 
may write S as a two dimensional curl of an integer func- 
tion n, S(r) = z X Vn(R), where n should be understood 

as residing on sites of the dual lattice, R (see Fig. 2). 
Thus we get an effective Hamiltonian for the /i — s- limit 



n{h ^ 0)vM = Wdg = ^ 5] |Vn(R)|2. 

R 



(6) 



Hug is the "discrete Gaussian model", which is one of 
the simplest solid-on-solid models of an interface^^. This 
is well-known to be dual to the two dimensional Coulomb 
gas, and so undergoes a Kosterlitz-Thouless transition; in 
the context of a solid-on-solid model, it is a roughening 
transition'^. It is helpful to recognize what is going on in 
the roughening transition in terms of the original S vari- 
ables before returning to the case h> 0. Since n(R) rep- 
resents the height of fluctuating columns in an interface, 
we recognize that bonds for which S(r) = z x Vn(R) ^ 
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actually represent domain walls between regions of differ- 
ent heights. Note that the condition V • S(r) = guaran- 
tees that the domain walls must be closed. The impor- 
tant point is that the roughening transition represents 
the temperature above which such closed domain walls 
proliferate, i.e., percolate through the system. When ex- 
pressed in terms of the Coulomb gas (i.e., vortex) model, 
the high temperature, rough phase corresponds to the 
low temperature, bound pair phase. 

Returning to the h > phase, we see that Eq. 5 al- 
lows us to work directly with the domain wall degrees 
of freedom. The fact that we allow configurations with 

V • S(r) 7^ means that we are allowing open as well 
as closed domain walls. In fact, locations for which 

V • S(r) ^ locate points where a domain wall comes 
to an end. This is precisely the situation one has when 
an interface contains screw dislocations: locations for 

V • S(r) ^ define the centers of these screw disloca- 
tions, and from the second term of TivM wc identify the 
core energy of the dislocation as l/2h. Remarkably, we 
have arrived at a theory that involves string-like objects 
(the domain walls) and vortex-like objects (the screw 
dislocations), which is very much in line with what we 
expected from the considerations discussed in the Intro- 
duction. We will see that Eq. 5 is essentially the dual of 
a model in which the string and vortex degrees of free- 
dom are explicit; however, we will need to introduce the 
vortex core energy before this is apparent. 

The above discussion shows that in our model (and its 
dual) the dislocations (vortices) are endpoints of domain 
walls (strings), so that when they are in a bound state, 
basic excitations of the system contain closed loops. We 
would like a representation of the model that captures 
this physics, presumably one that is closely related to 
Eq. 6. Toward this end we need to represent the domain 
wall degrees of freedom S as differences. This is most 
easily developed for the square lattice, so for the rest of 
this section and the next few, we specialize to this case; 
the triangular lattice will be dealt with in Section IV B. 
We start by writing 



5i(r) = mi{R = r + Ai/2 - A2/2) 

- mi(R = r + Ai/2 + A2/2), 
S2{r) = m2(R = r + Ai/2 + A2/2) 

- m2(R = r - Ai/2 -|- A2/2). 



(7) 



Note that we have located the integer fields mi, m2 on 
the dual lattice sites {R} of the original square lattice. In 
Sections IV A and IV B, we will see that representations 
such as Eqs. 7 are convenient for finding fixed points 

that respect the lattice symmetry. For now, we make 
the further transformation 

m2{X,Y)=n{X,Y), 

X-A2/2 

mi{X,Y)=n{X,Y)+ J2 M^'^^), (8) 

where R = {X,Y) are sites on the dual square lattice, 
n and A are integer fields, and the sum in the second of 
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FIG. 2. Real lattice sites (r) (black filled circles) and dual 
lattice sites (R) (open circles) of a square lattice. Ai = ooy 
and A2 = aox are the vectors denoting nearest neighbor 
bonds, and S2 arc bond integer fields that live on the 
nearest neighbor bonds as shown. Si and ^2 can be written 
in terms of integers mi and m2 on the dual sites {R} as de- 
fined in Eq. 7. A further transformation (Eq. 8) defines two 
integer fields n and A, where A lives on the vertical bonds 
and n resides on the dual lattice sites. 



Eqs. 8 is along rows in the square array, starting at some 
reference line xq which can conveniently be chosen as the 
boundary of the system. If one fixes the values of n and 
A along this boundary (e.g., n{xo,y) = 0, A{xo,Y) = 0) 
it is not hard to see that there is a one-to-one invertible 
mapping between the fields {81,82) and {n,A), so the 
new degrees of freedom capture all the domain wall con- 
figurations without overcounting them^". Notice in this 
representation it is natural to think of the n field as resid- 
ing on the dual lattice sites, while the A field resides on 
the vertical bonds of the original lattice. This is shown 
in Fig. 2. 

We now have the representation 



WvM = ^^|(Vn)r + A(a;,i/ + Ai/2)x| 

r 



dA\ 

dy). 



(9) 



In this representation, it is clear that the n field allows 
us to represent configurations with closed domain walls. 
The A field introduces the open domain wall configura- 
tions in two ways: it allows us to directly occ;upy the 
vertical bonds with non-zero values, and it allows us to 
remove the domain wall energy along vertical bonds of 
closed domain walls which allows open configurations 
with horizontal domain wall segments. Fig. 3 illustrates 
an example of this. 
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FIG. 3. A low energy configuration involving A 7^ 0. A re- 
gion of n = 1 (hatched squares) is embedded in a surrounding 
n = region (white squares) ; heavy line represents a domain 
wall. The line segment with A = 1 cancels the domain wall 
energy for part of its length, leaving an open domain wall. 



B. Vortex Representation 

To explicitly show the duality between dislocations and 
vortices, as well as to introduce the core energy of the 
vortices, we need a representation in which the vortex 
degrees of freedom are explicit. This is accomplished by 
eliminating n using the Poisson resummation formula, 
just as in the case h = 0. The partition sum becomes 



{"} {A} 
{A} {m}-^ 

- 27rz ^ 0(R)m(R) 



HvM[{n},{A}]J, 

-Wvm[W,{A}] 



R 



(10) 



The functional integral in Eq. 10 may be carried through 
if we represent the fields in terms of their Fourier trans- 
forms. Writing (j){q) = Er e"* "^i?i'(R-), with N the num- 
ber of lattice sites, and similar definitions for m(q), n(q), 
and ^(q), we have 



--VM 



{A} {m} 



q 



- |2 



^(/)(q)m(q) 



where Q = {Qx,Qy), and = 1 — e 



? Qy 



(11) 



1 



e "I'^i . Integrating out ^ gives, up to an unimportant 
prefactor. 



2yM 

= E E ^^p 

{A} {m} 



-E 

q 

N 
2m 



+ 



2K\Q\^ 2h 



E 



|w(q)P 



(12) 



N Sq ' 



Eq. 12 may be conveniently represented in real space 
in terms of Go{r) — ^ 
j_ 

N ■ 

that Go(r) 
now write 



'^I'VIQP and Gi(r) 

it is useful to nc 
logr and Gi(r) ~ arctan(y/a;). We can 



6 ^'^''^Qx/QijlQl'^- For large r, it is useful to note 



2vM 

= J2J2 ^^p| ~ ^TT^Jf J2 ™(R-)Go(R - R')w(R') 

R.R' 



{A} {m} 



dy 
OA 



2-EU^ G,(r-r>(R') 

r,R' ^ ^ ^ 



-E 

2h^ 



dA 
dy 



(13) 



The logarithmic interaction among the m variables estab- 
lishes the fact that they are the vortex degrees of free- 
dom. In this representation it is appropriate to introduce 
the core energy for the vortices, after which the resulting 
model is identical to the one analyzed in Refs. 1 and 8. 
The resulting system can be considered as a generalized 



Villain model, for which Zgvm = J2 



{m,m'} ' 



-Hgvm 



with 



Hgvm = 2Tr'^K ^ m(R)G'o(R - R')m(R') 

R,R' 

- 2m ^ m'(r)Gi(r - r')TO(R') 

r,R' 

+ ^E^'w'+^cE"^W- 



(14) 



In Eq. 14 we have replaced {dA/dy)r m'(r) to make 
the duality of the model apparent: upon interchanging 
Ec ^ l/2h and changing K l/4n'^K, the partition 
sum is unchanged. It immediately follows that if there is 
an unbinding transition for the vortices (the m variables), 
there must also be such a transition for the dislocations 
(the m' variables) at an appropriate location in the phase 
diagram. This duality was exploited in Ref. 8 to find 
the RG flow equations for large Ec and small h. In the 
limit analyzed below, large Ec and small l/2h, we shall 
see that there is a dislocation unbinding transition. The 
duality immediately tells us there must also be vortex 
unbinding for small Ec and large l/2/i. 
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While Hgvm displays an elegant symmetry, it is hard 
to work with directly. The appearance of an imagi- 
nary term in the Hamiltonian makes the configurational 
weights in the partition sum complex and can lead to sub- 
tle complications, particularly if one wishes to replace the 
integer degrees of freedom with continuous ones, which 
we will need to do to perform an RG analysis. Further- 
more, as discussed above the physics of vortex (disloca- 
tion) unbinding is profoundly affected by fluctuations of 
the strings (domain walls), which are now only implicit 
in Eq. 14. To see more clearly what is going on, it is 
worthwhile to return to a representation such as Eq. 9. 
We will proceed to do this, but in doing so we will also 
replace Hgvm by a continuum model that is amenable 
to an RG analysis. 



C. Continuum Model 

To return from the vortex representation to the domain 
wall-dislocation representation, one only needs to apply 
the Poisson resummation formula to the m sum in -Zgvm ■ 
However, for the RG analysis of this system we will want 
to develop an equivalent model that contains continuous 
rather than integer fields, and it is at this point that it 
is convenient to begin doing so. Whereas 



K 



eflF 



GVM 



= ^ ^ / expi -Hym {n <!>}, {m' ^ ^} 
-E^Y^n" (R) - 27ri ^ (/.(R)n(R) I , 

R R J 



{A}-' R n(R)=-oo I 

R R J 



dA 
dy 



is an exact representation of the generalized Villain 
model, we truncate the sums X^^r,)=_oo Sn(R)=-i- 
For large Ec this is an excellent approximation; vortices 
with large topological charge play little role in the prop- 
erties of the system when the fugacity is small. Using 
cxp(-27rm(/) - EcZIr"^) = 1 + ycos(27r^) w 
exp[j/cos(27rc6)] with y = 2er^'= twice the fugacity, we 
arrive at an intermediate model, 



-Z' = ^ y" expj - HvM [4>, ^ 



y y^cos(27rg 

R 



(15) 



At this point we have one continuous and one integer 
field. To go over to a fully continuous model we make 
the replacement Z' ^ Z = JV^JVa g-^efd'/'-a] , with 



/ 



-|V<A(r)+a(r)x| +-^- 



- y cos(27r0(r)) -|- y 
-yaC0s(27ra(r)) + 



(16) 



In Eq. 16, we have gone over from a lattice to a contin- 
uum representation, taking our lattice constant oq = 1 
as OUT unit of length, and we have subtracted constants 
so that the ground state energy is zero. Although our re- 
placement of the integer field A with the continuous field 
o is not a controlled approximation in the same sense as 
when we introduced (j) in favor of n, we have constructed 
Hcfi so that it has the same symmetry properties under 
translations of a ^ a + m, with m an integer, as had 
Wgvm, using the ya cosine term. Thus we expect HeS 
and TYgvm to have the same phases and types of transi- 
tions among them^. Our replacement of Hgvm with TieS 
is in fact no better or worse than our initial replacement 
of HxY with HvM- Ti-cS is the model that we will focus 
on in our RG analysis. 

It is useful to recognize that the low-energy configura- 
tions of {(f), a) in Ticff mirror those of the Villain model 
(Eq. 5) with which we started. For example, if we con- 
sider configurations with o = 0, then HeS is a two di- 
mensional sine-Gordon model, supporting kink excita- 
tions^^ that are directly analogous to domain walls. Hcs 
also supports solitons that represent screw dislocations. 
We can see this by considering configurations in which 
(j){x, 2/) — > 1 as X — > — oo and (p{x, y) — > as x — > oo. For 
a = this will force in a domain wall running along the 
y direction. However, if we set a{x,y) = 5{x), then we 
can produce a configuration of vanishing energy and sat- 
isfy the boundary condition on c6 with 4){x, y) = 0(— x), 
where G(— x) is a step function. A screw dislocation is 
forced into the system when we impose the boundary 
conditions [cf){x,y),a{x,y)\ ^ ['^SGs(a;), 0] as y ^ — oo, 
and [0(x, y), a(x, y)] [6(— x), (5(x)] as y ^ oo, where 
4>SGs{x) is the kink soliton of the sine-Gordon model^-'^. 
These boundary conditions guarantee that the domain 
wall must end somewhere in the bulk of the sample. This 
endpoint is the screw dislocation core. This is pictorially 
shown in Fig. 4. 

Finally, while most parameters entering HesiK, h, and 
y) are directly related to parameters of the generalized 
Villain model {K, h, and Ec)-, the value of ya in terms of 
these parameters is not immediately obvious. The con- 
nection can be made by recognizing that ya plays a key 
role in determining the dislocation core energy in TieS- if 
ya is excessively large, then variations of a in the y di- 
rection become energetically very expensive, so that the 
core of a dislocation, where a for example varies from 
to S{x), becomes high in energy. 

Unfortunately the screw dislocation is a rather com- 
plicated soliton so that it is not obvious how one might 
obtain an analytic expression for its energy. We can, 
however, make a rough estimate for ya by looking at the 
dislocation configuration in the Hamiltonian implicit in 
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(a) Closed wall configuration: 




(P(x) 



Sine-Gordon soliton 



(b) Open wall configuration: 



a =0 



,-i^|a=5(x)| / 










a(x)=6(x) 


9 ' 





a=0 



a=0 



(p(x) 



a=0 



FIG. 4. Difltoront low energy soliton configurations of {4>, a) 
for the effective Hamiltonian TieS defined in Eq. 16. (a) 
For closed wall configurations with a{x) = 0, TieS is a 
two-dimensional sine-Gordon equation with sine-Gordon soli- 
tons, (b) a{x) = 5{x) for large positive y, a{x) = for large 
negative y, and 4>{x,y) —^laax—^ —oo, (j){x,y) ^ as 
a; — > oo. The domain wall has finite energy per unit length 
as ?/ — > — oo, but vanishing energy per unit length as j/ — » oo. 
The resulting configuration is a screw dislocation. 



the partition sum for Z' in Eq. 15. In Eq. 15, A is an inte- 
ger variable residing on the vertical bonds, and we create 
a dislocation near the origin by setting A{x,Y) = —1 for 
a; = 0, Y > 0, and A = elsewhere, and introducing the 
same boundary conditions on ^ described two paragraphs 
above. It is clear for large negative Y, cj) will contain a 
sine-Gordon kink soliton, whereas for large positive Y, 
it has the form of a step function. The (j) kink presum- 
ably narrows from the sine-Gordon width to zero in the 
vicinity of the origin. The energy of the configuration 
will take the form e = + l/2h, the first contribution 
representing the energy stored in the </> field and its in- 
teraction with A, and the second the energy cost for the 
step down in A at the origin. In going over to Ti.cS, where 
A is replaced by the continuous field a, if we assume the 
contribution is unchanged, then we can approximately 
match the dislocation core energies by matching the cost 
of the step down in a with l/2h. This is estimated by 
keeping only the second and fifth terms of the right hand 
side of Eq. 16, so that the step is again a sine-Gordon 
kink. The energy of the kink is^^ 8y^ya/h; we obtain the 
match by setting ya ~ l/256/i. 

With this estimate, it is possible to connect the phases 
of HeS with those of the GVM and XY models. We now 
address this with an RG analysis. 



-71 



I 



.71 



Adl 



FIG. 5. Momentum shell in reciprocal space for the square 
lattice. Ai and A2 are equal length vectors of magnitude tt 
(ao = 1) to the center of the BZ edges. The width of the 
shaded region is AM. 



III. RENORMALIZATION GROUP ANALYSIS 

A. Derivation of Scaling Relations 

Our analysis of HeS proceeds with a momentum shell 
renormalization group procedure'^ that is perturbative in 
y and ya- This involves dividing the fields into long and 
short wavelength components: ^(r) = 0^(r) -|- <^^(r), 
with 



M,\qy\<A/b 



(27r)2 



</.>(r) 



/ ) 

JA/b<\qa,\<A JA/b<\qy\<Aj 



e-^^-'^^iq). (17) 



There is an analogous decomposition for a. In Eqs. 17, 

A = Tv/ao, b = is the rescaling factor, and we are 
including in <p^ momentum components in a square shell 
at the edge of the Brillouin zone, as illustrated in Fig. 5. 

As we shall sec, the fixed points we are interested in are 
not orientationally isotropic, so it is important to choose 
the momentum shell so that it respects the symmetry 
of the lattice. (When we address the triangular lattice 
below, our momcntmn shell will be a hexagon.) 

As usual the analysis proceeds by integrating out (j)-' 
and from the partition function Z. For y, ya = we 
then have 



Z{£) = e-^W J Ty^< j 

l_ f< d^q 
KJ 



X exp 



2K J (27r)2 
+ e\a<{q)? 



-«q^ (q) +a<(q)x 



- |2 



(18) 
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where e is a prefactor that comes from the integrals 
over (jr' and a^, /"^ indicates an integral in the inner 
part of the Brillouin zone (i.e., inside the shaded region 
in Fig. 5) and ^ — \/ K/h is the width of a string con- 
necting a highly separated vortex- antivortex pair-'^". We 
can bring what remains of the Hamiltonian back to its 
original form (for y, ya — 0) with the scaling transfor- 
mation 



q = q'/ft , r = br' , 
i<(r) = (A'(r') , a<(r) = a'(r')/6. 



(19) 



We then have in real space 

Z(£) = e-^W Jv<j)'Jva' 



X exp< - — / d'r' 



+ e\a'{r')\^ 



|V</.'(r')+a'(r')x| 



with K' = K, ^'^ = For a very narrow shell 

(£ ^ di, 6 — > 1 + di), the differential changes in K and 
are just 



dK 

Hi 

de 

de 



0, 



(20) 



Thus we have a set of fixed points whose stability with 
respect to the cosine perturbations we need to assess. For 
the y cos{2'K(f)) perturbation we will proceed in the stan- 
dard fashion"^'"*, and we will see near the fixed points of 
interest it has a large anomalous dimension - the pertur- 
bation is highly irrelevant. The ?/cos(27ra) perturbation, 
on the other hand, is rather unusual. The scaling rela- 
tion for a in Eq. 19 causes it to shrink as one integrates 
out the short wavelength fluctuations of the system. This 
suggests we should expand 



yaC0s{2na) - ya 



^ 2/2)1 



- (2n)! 



(-1)"(2^«) 



2n 



(21) 



with the initial values y2n{l = 0) = Va- Our effective 
Hamiltonian is now 



Weff = / d'r 



-|V0(r)+a(r)x| + - 
y cos(27r0(r)) -|- y 



y2n 



^ (2n)! 

n=l ^ ' 



(-l)"(27ra)2 



(22) 



If we simply substitute a(r) — > a'{r')/b in the last term 

of Eq. 22, we sec that the coefficients y2n ^ This 
suggests, as we will confirm below, that most of the terms 
in this expansion are irrelevant. However, the Gaussian 



(i.e., first) term is special because it neither grows nor 
shrinks, suggesting we should incorporate it into our fixed 
point. Writing p{£) = 47r^y2(^), we take our unperturbed 
Hamiltonian to be 



Ho 



-I 



da 
dy 



a(r)x|^ 



(23) 



To first order in y,ya, the partition function takes the 
form 



>e-Wo 



Z = J V(l)< j D(/>> j Va< j Va 

X ll-y j dVcos[27r(?!)<(r) + 27r(/)>(r)] 



- / d'r 



OO 

E 



y2n 



(-1)" 27ra<(r) -h27ra>(r) 



(24) 



The "1" term above has essentially been discussed al- 
ready; the only difference one gets after integrating out 
(f)^, a^ from Eq. 18 is that F{£) now has a correction of 
order ya due to the p term appearing in Hq. Thus the 
scaling relations Eqs. 20 are correct to first order in y 
and ya- We next wish to integrate out the short wave- 
length fields in the middle term of Eq. 24. Calling this 
contribution Zy, the integral is 



-y j V<f)< j V<j)> j Va< j Va 
rf^rcos[27r0<(r) -|- 27r(;&>(r)] 

-F(t) 



>e-«o 




4/ 



d\ 



1^- 



l + Kp + eql 
J V4>< J Va<e-'^o J (frcos\2Tr(f)<{r^ 



(25) 



where Ti.Q is obtained by Fourier transforming the right 
hand side of Eq. 23, and then dropping wavevectors in the 

resulting q integral that are in the momentum shell. We 
can already see the crucial role played by the parameter 
p: for p = 0, the d^q integral is divergent, and it is 
apparent that Zy vanishes in the thermodynamic limit. 
Since the values of p in which we are interested are indeed 
small, it is apparent that the j/cos(27r0) is going to be 
strongly irrelevant. We estimate the / d^q integral by 
evaluating it for ^ = since this variable is irrelevant 
(Eq. 20). For small p, one finds 



r 



d\ 



l + Kp 



= ttK 




\ogb + 0{p). 



Rescaling the fields and lengths, we can return Zy to its 
original form, writing 



Zy = -y'e 



-Fit) 



j V(j) j Va e"^" j rf^rcos[27r(^(r)]. 
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with 



y =yb W '<p , 



where again we arc focusing on small p. Recalling b = 
1 + 0{di'^), we now have the scaling relation for y, 



de 



ttK. 



'l + Kp 



Kp 



(26) 



Because p is small (recall its initial value is 0{ya), and 
we shall see it shrinks under renormalization), for the 
parameters of interest y shrinks very rapidly as we in- 
tegrate out short wavelengths. The ycos{2iT(j)) term is 
thus strongly irrelevant, suggesting that the system is in 
a "rough" phase - domain walls are proliferated through 
the system. We will discuss the implications of this be- 
low. 

We are left finally with deriving the scaling relations 

for the polynomial terms in a. These are also handled 
in a standard fashion"'. We first consider a term in the 
last sum of Eq. 24. Expressed in terms of the Fourier 
components of the fields, such a term may be written as 



(2n)!^"^^"/ d\id\2...d'q2. 



2/2. 



X (27r)^J(qi + q2 + ■ . . + q2n) 
X a(qi)a(q2)...a(q2n)- 



(27) 



In writing this in terms of short and long-wavelength con- 
tributions, we wish to separate each of the a(qi)'s in the 
product above by the location of qi in the Brillouin zone. 
One thus has 

a(qi)a(q2) • • •a(q2n) 



+ 



a<(qi)a<(q2) . . .a<(q2 
2n(2n- 1)^> > 



a>(qi)a>(q2)a<(q3) • . .a<(q2„) + (28) 




-I- 



FIG. 6. Pictorial representation of Eq. 29. Integral over 
short wavelength fields [4>^ (r)] of the y2n term introduces a 
renormalization of the j/2n-2 term. The graphical decompo- 
sition of a 2n point vertex is shown. The first diagram on 
the right hand side is a tree diagram that represents all the 
long-wavelength fields [^^(r)] which correspond to the first 
term in Eq. 29. Contraction between a pair of legs is shown 
in the second diagram and it represents the second term in 
Eq. 29. The black dots in the figure represent other legs. 



'^''^ a>(-q)a>(q). (30) 



(2^)2 



Eq. 29 shows that the integral over short wavelengths of 
the 2/271 term introduces a renormalization of the y2n-2 
term. This is illustrated graphically in Fig. 6. Rescaling 
fields and lengths, we see that Zy^^ may also be returned 
to its original form, 



-F(l) 



Va e 



-Ho 



-1)" 



(2n)! 



27ra(r) 



2n 



(31) 



with 



V2n 



y2nb 



2-2n 



27r^[£{p,0d£] y2. 



:ri+2- 



This allows us to write down our last scaling relation. 



dy: 



2n 



dl 



= -(2n - 2)y2n - 2t,'^ C{p,0y2n+2. 



(32) 



A slightly subtle but important question arises when we 
where wc have dropped the term with just one a>, an- recall 47r^y2 = P- Should Eq. 32 be applied to p, even 



ticipating this will integrate to zero in the functional in- 
tegral. Because our momentum shell has width d£, we 
need only retain the two terms explicitly shown in Eq. 28; 
terms with more factors of involve higher powers of 
d£ and vanish when we take d(. — > 0. Substituting this 
expansion into Eq. 27, we obtain 



Zy^n = -y2n I V<p^ I Va<e-'^o 



X e 



though it appears in Tig and is not included in the Zy^^ 's? 
The answer is yes; the integration of the term does 
produce a Gaussian term that renormalizes the ^pa^ in 
Ho . This is easily seen if we reexponentiate the 0{ya) 
contribution to the partition function. The renormaliza- 
tions of the y2n terms follow through exactly as wc have 
described above, and one can then see explicitly that the 
1/4 term renormalizes p. Thus we can write for the n = 1 
case of Eq. 32 



a, (-1)" 



d'r 



(2n) 



27ra<(r) 



2n 



(—l)n-l r- -|2ri-2| 

^^-L-^27r'jC{p,0di[2na<{v)\ \, (29) 



dp 
~dl 



-^ = -87r^£(p,C)2/4 



(33) 



with 



The derivation of the scaling relations is completed by 
writing down an explicit expression for C{p, ^). The func- 
tional integrals in Eq. 30 yield 
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f 



2Kq^ 



(34) 



This integral can be computed exactly in a tedious but 
in principle straightforward calculation. The result is 
displayed in Appendix A. For small p, the result (Eq. 60) 
may be expanded to yield 



(35) 



Eq. 35 once again demonstrates the importance of the 
parameter p: as p ^ Q, L diverges, suggesting that the 
\pa^ term in Ho is irrelevant if the initial value of p is 
small enough. We will see this is indeed the case. 

Eqs. 20, 26, 32, and 33 arc the scaling equations for 
our problem. Before presenting results obtained from 
numerically integrating them in Sec. V, we discuss the 
fixed points of the equations and their interpretation. 



B. Fixed Points: Bound and Unbound Dislocation 
States 

It is easy to see that fixed points of Eqs. 20, 26, 32, 
and 33 occur for y = 0, y2n>i = 0, and = 0, so that 
the generic fixed point Hamiltonian has the form 



1 
2K 



|V0(r) 



■ a(r)x|^ 



1 



pa(r)^ 



• (36) 



We will see in our numerical integration of the scaling 
relations that any (small) value of p is an attractive fixed 
point. This will give us an unbroken line of fixed points. 
In spite of the p = point being connected to the p > 
fixed points, these two cases are physically different. The 
reason goes back to the role a plays in allowing us to form 
separated dislocation pairs. We demonstrated in the last 
section that one can remove the energy of a section of 
domain wall in the (p field when p = 0, by creating a 
narrow line segment where a ^ along some part of the 
domain wall (see Fig. 4). Because the t/cos(27r0) term is 
irrelevant, we understand that domain walls in cj) prolifer- 
ate throTigh the system. Thus, if we consider only a = 
configurations, we can find arbitrarily large domain wall 
loops. Such loops can be broken open to create disloca- 
tion pairs. The only extra energy cost for doing this to a 
closed loop of length L is twice the dislocation core en- 
ergy; the loop can be deformed to move the dislocations 
a distance ~ L apart, without any further cost in energy. 
This finite core energy is more than balanced by the en- 
tropy associated with the number of open domain wall 
configurations for large L, so that large open domain wall 
configurations - highly separated dislocation states - are 
not suppressed by the dislocation core energy. We thus 
identify the p = fixed point as an unbound dislocation 
state. 

For any non- vanishing value of p, the situation for large 
enough L becomes profoundly different. The "mass" 



term for a tells us that we cannot completely eliminate 
the domain wall energy. If we create a dislocation pair 
by deforming a closed loop, there is a residual attrac- 
tive interaction between them. We can thus identify p 
as a renormalized string tension for the dislocation pair. 
When p > 0, the dislocations are bound into pairs. It is 
interesting to note that because the a field lives only on 
the vertical bonds of the underlying lattice, we do not 
get the configurational entropy that the domain walls in 
the (f) field has - so even though the residual "string" 
connecting the dislocations may have much lower energy 
per unit length than the closed domain wall loops, the 
free energy of closed domain walls can vanish, while that 
of the residual strings remains finite. 

While it is tempting to interpret the energy of the 
residual strings as leading to linear confinement, this is 
not correct because it ignores the effect of fluctuations 
in (p. A simple interpretation may be had if we interpret 
p > as telling us that the"integerness" of the original 
A field cannot be ignored - the cos(27ra) term contains 
physics which is relevant at long wavelengths - even while 
the irrelevance of y tells us we may treat (j) as a contin- 
uous field and ignore the cos(2TT(j)) perturbation. This 
suggests we can return to Z' (Eq. 15), setting y = 0, 
which allows us to integrate out (j). Up to unimportant 
prefactors, one finds Z' oc ^^a} e"^*^*^, with 



dA 
dy 



Recalling that {dA/dy) is a representation of the dislo- 
cation field, we see that we have come to a Coulomb gas 
Hamiltonian! This suggests we should think of the p > 
state as a set of logarithmically interacting dislocations, 
rather than linearly confined ones. However, in coming 
to HcG wc understood p > as signaling that we cannot 
ignore the discreteness of the underlying charge. This is 
the usual situation when the charges are in a bound state. 
If one can treat (dA/dy) (da/dy) as a continuous field, 
the state may be interpreted as a "charged liquid" and 
the dislocations may be treated as unbound. This is only 
valid if p = 0. 

An alternative look at the fixed points may be had 
by reexpressing them in terms of the vortices. To do 
this, we return (p to its integer form n but keep a as 
a continuous field, adopting the fixed point form for its 
energetics. Applying the Poisson resummation formula 
to n brings us to an effective partition function Z' = 
Ew/2^a e-«' with 



1 



2K\Q\ 



2Tr''K^^ 27ri|;^a(q)TO(-q) ^. 
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Integrating out a gives, again up to unimportant prefac- 



tors, Z' = J2 
now 



^, with the reduced Hamiltonian H 



n 



q 



N 



\Q.\ 



\Qy\^+pK\Q\^ 



The energetics of the vortices is profoundly different for 
p > and /9 = 0. This is most clearly seen by considering 
a configuration for a single vortex at the origin, m(R) = 

<5r,0 5 in a finite size system of linear dimension L. The 
energy of such a configuration has the form 



<H> 



Iv 



(1 



C)logL, 
iP = 0) 



(p>0) 



(37) 



where C = j^"" dO [cos^ 6l/(sin2 6l + pK)]. We find the 
usual logarithmic energy for a vortex when p > 0, al- 
though the correction C to the coupling constant be- 
comes very large for small p. The effect of the dislo- 
cations thus is a strong upward rcnormalization of the 
effective K in this state. For p = 0, by contrast, the en- 
ergy of a single vortex configuration grows linearly with 
L. Thus, the state that we identified as the unbound 
dislocation state becomes, in the dual representation, a 
linearly confined vortex state. By contrast, the p > 
state represents logarithmically confined vortices, or, in 
the dual language, logarithmically bound dislocations. 
Apparently these dual defects are in the same state at 
such fixed points. 

The duality of the partition function tells us then 
that there are three possible states for the vortices, lin- 
early confined, logarithmically bound, or deconfined - 
the last because the symmetry of the effective Hamilto- 
nian Hgvm (Eq. 14) tells us that if a deconfined dis- 
location state exists, so must a deconfined vortex state. 
Recent simulation studies^^ have strongly supported the 
existence of three such states in the XY model. For 
example, if one fixes K below the critical value for the 
Kosterlitz-Thouless transition, then at small h one finds 
fluctuations that scale with system size i in a way ex- 
pected for unbound vortices. At intermediate values of 
h, the fluctuations behave as if the vortices are in pairs 
that are interacting logarithmically at large separations. 
For larger values of h, they behave like dipoles with linear 
confinement^^. 

Finally, one may notice that H is not invariant under 
Qx Qy, so that the fixed points we are examining do 
not retain the full symmetry of the underlying lattice. 
Indeed, the loss of symmetry was been apparent in our 
development of the RG scaling relations, for which we 
saw Tio and ultimately has configurations of vanish- 
ing energy along the qy — axis, but not on the Qx ^ 
axis. This problem originated when we went from Hgvm 
to our HeS- For the former model, the energetics of a dis- 
location with a domain wall exiting along a vertical bond 
is the same for one in which it exits along a horizontal 
bond; for the latter, they are different. The model can 
be formulated in a way that retains the lattice symmetry. 



at the sacrifice of some of the simplicity of the model we 
are currently analyzing. We will detail how this is done 
below. For now we note these models generate unbound 
dislocation (and vortex) fixed points that have not one 
but several nodal lines in the Brillouin zone, in a way 
that respects the orientational symmetry of the lattice - 
so that the H* we are currently analyzing is really only 
one of a class of fixed point Hamiltonians that can arise 
in XY spin models with symmetry-breaking fields. 



C. Integration of Scaling Relations 

We now need to show that there are parameter regimes 
for which TicS flows to both with p = and p > 
0. The first two scaling relations - Eqs. 20 - are trivial 
to deal with. The first of these simply states that K 
is invariant as we integrate out short wavelengths: we 
simply take this to be a constant in the discussion that 
follows. The second is easily integrated to give 



-It 



(38) 



where = \/K/h is the bare string width. Because of 
the simple relation between ^ and i, it is convenient to 
characterize the length scale to which we have integrated 
in terms of ^; i.e., b = e^ = £,o Thus we will present 
RG flows below in terms of ^. 

Eq. 26 is also easy to deal with since the parameter y 
docs not rcnormalize any of the ?/2ri's. Since the initial 
value of p is small in our perturbative approach, and it 
shrinks under the RG, y{t} shrinks rapidly as ^(f) 0. 
We will not show this explicitly. 

Eqs. 32 by contrast appear formidable since they 
present an inflnite tower of equations. However, because 
the shell integral C is the same for each of these, the 
solutions of the equations are related to one another by 
y2n+2 = e"^^t/2n = ?^2/2n/'?0' ^s may easily be conflrmed 
by direct substitution. Recalling p = 47r^t/2, this tells us 
Eq. 33 may be written as 



dp _ _2TfXp^ 
d£- 

Finally, noting d^^/di = -2^2^-2^ 
rewrite this as 



dp 
d^ 



'^^C{p,Op 



-2^2(f), we can 



(39) 



Eq. 39 is particularly convenient for integration. Note the 
explicit dependence on shows that the rcnormalization 
of p is strongest for the largest values oih: p scales to zero 
for large values of h, leading to the unconfined dislocation 
phase. 

Fig. 7 shows a typical result for the integration of 
Eq. 39. For these flows, we have used our estimate of 
ya in terms of h to fix the initial value of p, so that all 
the flows begin from the same point in the plot. The 
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FIG. 7. RG flows for Eq. 39. Left vertical axis is a fixed 

line, and heavy line is a separatrix between flows that reach 
p = for flnite / and flow to the origin as ? ^ oo, and those 
that have p > at the end of their flow. Finite values of p in 
Ht, indicate dislocation binding. 



structure that was described in the Introduction is ap- 
parent: there is a fixed line of p values, with the p = 
point - the unbound dislocation phase - an endpoint of 
this line. For small values of ^o, p(^) lies below the sepa- 
ratrix (shown as a heavy line), touches the p = axis at 
a finite value of ^^/^q = e~^^ , and remains on this axis 
as it flows to the origin. Note that the divergence in C as 
/) ^ guarantees the flows never cross the /? = plane 
to negative values, which would indicate an instability in 
our model. We can interpret £scr = ao^^ as a screening 
length: for separations below this the dislocations appear 
to be bound, while for larger separations the interaction 
is screened by other dislocations, allowing for an unbind- 
ing transition. Apparently iscr diverges as the transition 
is approached. This means the deconfinement transition 
has a feature in common with other more standard phase 
transitions: there is a diverging length scale. Above the 
separatrix, the flows end at a non-vanishing value of p, 
which as disciisscd above indicates the dislocations are 
in a bound state. The separatrix represents a deconfine- 
ment line, with dislocations unbound if is smaller than 
some critical value for a fixed valuc^ of K. 

The numerical integration of the scaling relations allow 
us to find some critical exponents for how p develops and 
how £scr diverges as the phase boundary is approached. 
Before giving results for these, we will demonstrate that 
the lattice symmetry can be properly incorporated into 
the model and find the scaling relations analogous to 
those above. This is the subject of the next section. 



IV. FORMULATION OF THE MODEL WITH 
FULL LATTICE SYMMETRY 

A. Square Lattice 

In this section we will formulate the theory retaining 
the full underlying symmetry of the square lattice. We 
begin by defining the bond variables 5*1 and 5*2 in terms of 
integer fields mi and TO2 on the dual lattice sites {R} as 
in Eq. 7 and substitute into Eq. 5 to obtain in momentum 
space 



q 



2 I 1 2' 

+ Q2TO2(q) 



+ ;^5Z^I<5l<32ml(q)+QlQ2m2(q)|^ (40) 
q 

where Qi = 1 - e-^i'^^ and Q2 = 1 - e-'^J'^^ As 
previously, to carry out an RG analysis we need to replace 
'iTii^2 by continuous fields (^1,2 and add in cosine terms 
which tend to preserve the "integerness" of the fields, as 
well their sums and differences. Most such terms work 
out to be irrelevant; however, one particular combination 
turns out to be analog of the j/a cos(27ra) term of the last 
section. This is 



l-cos|27r [(^i(r + A2)-</)i(r)] 



+ 1 - cos<^ 27r [(/.2(r + Ai) - ,^2(r)] 



(41) 



It is useful now to do a Fourier transform of the integer 
fields mi and m2 in Eq. 41 and expand the result in the 
small q limit 



^i(r + A2) - 01 (r) ^ 1 ^(iq. A2)e''»--</.i(q), 

q 

+ Ai) - Mr) « ^ • ^iV^-'M^)- 

q 

The effective Hamiltonian is now 

^eff = ^ • ^^f\M<l)\' + (q • A2)^|<^2(q)|^ 

q I 

+ ^|(q- Ai)(q- A2)0i(q) 
+ (q- A2)(q- Ai).^2(q) |'| 

2tt 



cos 



N 



^(iq- A2)e''»--<^i(q) 



-|- I 1 — cos 



-5^(iq.Ai)e^^-</)2(q) 
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+ 1 - cos[27r^!>2(r) 



(42) 



In coming to Eq. 42, wc have taken the long wavelength 
limit Qi = zq • Ai and Q2 = iq. - A2. We note the sym- 
metry of Tigff in Eq. 42 under x y S'Hd 'Pi </'2 
preserves the square symmetry. 

Wc arc now set to perform the momentum shell renor- 
malization group procedure on the above effective Hamil- 
tonian that is perturbative in tja- Note that the y term is 
irrelevant and shrinks to zero as discussed in Sec. IIIC. 
Our analysis that follows is similar to that described in 
Section III - we divide the fields into long and short 
wavelength components (l)f and (f>f, integrate out the 
short wavelength fields (within the momentum shell as 
shown in Fig. 5) which give an overall prefactor in the 
partition function, and finally rescale the fields as in 
Eq. 19. In doing this we expand the ya cosine terms in 
Eq. 41 and incorporate their Gaussian contribution into 
our unperturbed Hamiltonian, which now becomes 



+ ±{{ci.A^f+pK{ci.A,f}\Mci)\' 



2K 



(q- Ai)(q- A2)0i(q) 

2 



+ (q- A2)(q- Ai)02(q) 



(43) 



where p{£) = ^TT^y2{i) and ^{i = 0) = ^/KJh. After 
rearrangement, the resulting unperturbed Hamiltonian 
takes the form, 





(q-Ai)^ 


+ pK{q-A2f 




Ai)^(q 


A2)^}|<^i(q)|^ 




(q-A2)2 
A2)'(q 


+ pK{ci-A^f 
Ai)'}|<^2(q)|' 



+ ^(q- Ai)2(q. A2)"<^i(q)</.2(q) 



(44) 



As in Sec. Ill A, the fixed points of the RG equations 

[Eqs. 20, 26, 32, and 33] occur for ?y2n>4 = 0, and = 0, 
so that the generic fixed point Hamiltonian is 



— {(q.Ai)2+pi^(q.A2)2}|</.i(q)|^ 

(45) 



— {(q.A2)^+pi^(q-Ai)2}|,^2(q)| 



One may notice that the fixed point Hamiltonian in 
Eq. 45 is invariant under Qy which implies that the 



fixed points under investigation do retain the full under- 
lying symmetry of the lattice. This is different than the 
case we studied in Sec. IIIB. As before, the scaling rela- 
tions show that any small value of p is an attractive fixed 
point, we identify the p = fixed point as the unbound 
dislocation state, and p > as logarithmically bound dis- 
location states. Because of the full lattice symmetry of 
TC* there arc now two nodal lines in the Brillouin zone 
when p — 0. 

For our later calculations we will ignore the cross terms 

in Ho ((/>i</>2) as those terms contribute only at higher 
orders of (starting from ^^), and the calculation is 
greatly simplified because {(j)i(t>2)o = for the Hamilto- 
nian in Eq. 44 when this is done. With this simplification, 
the propagators for the fields (f)i and ^2 calculated from 
Eq. 44 are 



1 



(q • Ai)2 + pKiq ■ A2)2 + ^^(q • Ai)2(q • A2)2 



^2(-q)02(q))o 



(q • A2)2 + pKid ■ Ai)2 + e2(q . Ai)2(q • A^)^ 



■ (46) 



With these we define the momentum shell integrals 
Cq {P, ^) for the the square lattice 

Ch{p,^)de = K y'^^(q. A2)2(</.i(-q)<^i(q)), 



Cl{p,OM = K 



(2^) 



(27r)^ 



(q- Ai)2(02(-q)02(q)). (47) 



Note £□=£□= Cu because of the square symmetry. 
The scaling relation in Eq. 39 remains the same, provided 
we substitute £ ^ £□ + £□ = 2Cn- The shell integral 
Cn may be computed to yield 



-Cd(p,6 



+ 



pK + 

[pK + ^2A2)3/2 



: arctan 



7r^l+eA2 



arctan I ^ J . (48) 



where we have returned to unitless distances, Ai_2 = 
ao = 1. 

For small values of pK we perform a Taylor series to 
obtain 



2n^pK{l + e2A2) ■ 



(49) 



We thus sec that the scaling relation for p in this for- 
mulation is identical to Eq. 39. The only caveat is that 
in our previous formulation, this relation held to all or- 
ders in £2^ whereas here there are corrections of C'(£^) 
and higher. In Sec. V we will present further results that 
follow from numerically integrating the scaling relations. 
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FIG. 8. Real lattice sites (r) (black filled circles) and bond 
center sites (patterned circles) used to represent Si , S2 and 
S3 as differences for a triangular lattice. Ai, A2 and A3 arc 
vectors denoting nearest neighbor bonds. Note that the lines 
connecting the m's cross the bonds they define, but they are 
not perpendicular to them. 



B. Triangular Lattice 

Wc now consider our theoretical model on a triangular 
lattice, following the same procedure as in Sec. IV A. 
For our triangular lattice we have three bond variables 

Si, 82, and 5*3 which are defined via three integer fields 
mi, 7712 and TO3 (respecting the symmetry of the lattice) 



(50) 



5'i(r) = mi(r + A3/2 + Ai) - mi(r + A3/2), 
S2{r) = m2(r - Ai/2 + A2) - m2(r - Ai/2), 

Sair) = m3(r + A2/2) - m3(r + A2/2 - A3). 

One may note that in Eqs. 50 the integer fields mi 2,3 
are not defined on the dual lattice sites as was done for 
the square lattice: in this case they are on the centers of 
the bonds, as illustrated in Fig. 8. The lines connecting 
the m's cross the bonds they define, but they are not 
perpendicular to them. The direction of the Ai's are 
also shown in Fig. 8. Substituting Eqs. 50 in Eq. 5 we 
obtain the unperturbed Hamiltonian in momentum space 

HyM 

= Si + |<92m2(q)|^ + |(33m3(q)|^ 

1 , ,2] 
+ ^|QiQ2mi(q) +(32(33m2(q) +<33<3i»n3(q)| >, (51) 

where Qj = 1 — e~*'i '^j. As before we replace the in- 
teger fields mi_2,3 with continuous fields (/'i,2,3 and add 
in cosine terms. The cosine term that is analogous to 
j/aCOs(27ra) is 



^1-C0s|^ 



W = > ^ I 1 - cos|27r [<^i(r + A2) - <^i(r) 
1 - cos<l 27r [02 (r + A3) - (j)2{r) 



+ I - cos<^ 2tt [(j)3{r + Ai) - <j>3{r 




(52) 



We now do a Fourier transform of the fields mi, 2. 3 and 
expand the result in long wavelength limit 0i(r + Aj) — 
Mr-) ~ ^Eq(^q- Aj)e^''--<^i(q) where i j = 1,2,3. 
The resulting effective Hamiltonian is 



N ^ \ 2K 



(q-Ai)2|</,i(q)|2 



+ (q-A2)^|<^2(q)|' + (q-A3)2|<^3(q)|' 



+ ^|(q-Ai)(q-A2)(^i(q) 

+ (q • A2)(q • A3)<^2(q) + (q • A3)(q • Ai)03(q) 



cos 



-5^(iq.A2)e^^-,^i(q) 



+ I I — cos 



+ I I — cos 



N 



^^(iq.A3)e^^-02(q) 



N 



27r. 



^■^(iq.AOe^^-03(q) 
q 

,^|(^I-cos[27r.^i(r)]^ 
+ { I - cos[27r(/)2(r)] ) + ( I - cos[27r(/)3(r 



In deriving Eq. 53 we took the long wavelength limit 

Qj = zq • Aj (j = 1, 2, 3). We may note that the above 
effective Hamiltonian is symmetric under interchanges of 
the (1,2,3) indices, thus preserving the triangular lattice 
symmetry. 

Our next step is to do the momentum shell RG on 
the effective Hamiltonian. As the procedure is identical 
to that discussed in Section IV A we present the main 
results. Expanding the cosine terms and absorbing the 
quadratic term into the unperturbed Hamiltonian, one 
finds 

Ho = 



-E - 



N 



2K 



{(q-Ai)' + (pif)(q-A2)2}|</.i(q)|2 

+ {(q-A2)2 + (pi^)(q-A3)2}|02(q)P 
+ {(q- A3)2 + (pJr)(q- Ai)'}|03(q) 
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+ 



2K 



(q- Ai)(q- A2)(/)i(q) 
+ (q- A2)(q- A3)02(q) 

+ (q- A3)(q- Ai)^3(q) 



Rearranging the terms in Eq. 54 we get, 

q I L 

+ e2(q. Ai)2(q. A2r}|^i(q)P 

+ {(q- A2)'+pif(q- As)^ 
+ ^2(q.A2)2(q.A3)2}|</)2(q)|' 
+ {(q-A3)^ + p/^(q-Ai)' 

+ e'(q-A3f(q-Aif}|<^3(q)|^ 



(54) 



K 



(q • Ai)(q • A2)2(q • A3)0i(q)02(q) 
+ (q • Ai)(q • A2)(q • A3)202(q)<^3(q) 
+ (q • Ai)2(q . A2)(q • A3)<^i(q)</>3(q) 



(55) 



Once again, the generic fixed point Hamiltonian is 

2/2n>4 = 0, ^2 = 0, and 



N ^2K 
q 



{(q-Ai)2+pi^(q.A2)2}|,^i(q)|2 

+ {(q-A2)2+pK(q.A3)2}|02(q)P 
+ {(q-A3)^+pi^(q-Ai)'}|<^3(q)|' 



(56) 



Because the fixed point Hamiltonian in Eq. 56 captures 
the triangular symmetry of the lattice, there are three 
nodal lines in the Brillouin zone for p = 0. As before 
we identify this as the unbound dislocation phase, p > 
again corresponds to a bound dislocation phase. 

For simplicity of our further calculations we will ignore 
the cross terms {4>i4>2T4>24>i,<t>i4>i) in as those terms 
contribute only at higher orders of In that case the 
propagators of Eq. 55 can be easily evaluated, yielding 



(0i(-q)0i(q))o 



(q • Ai)2 + pK{^ ■ A2)2 + C2(q • Ai)2(q • Aa)^ ' 



((/)2(-q)02(q)>o 




FIG. 9. Momentum shell of width KM in reciprocal space 
for the triangular lattice. Ai,A2, and A3 are equal length 
vectors of length 27r/(V3). Ai, A2 and A3 represent direc- 
tions of bond vectors of length ao{= 1) in Fig. 8. 



(</>3(-q)</>3(q))o 
_ 1 

~ (q • A3)2 + pi^(q • Ai)2 + ^2(q . A3)2(q . Ai)2 " 

The momentum shell integrals that we need to evaluate 
for the scaling equation are 

C\M = K |^^(q.A2)^</.i(-q)</.i(q)), 

£id^ = i^y"'^(q.A3)^(.^2(-q).^2(q)), (58) 

> d2q 



C\dl = K 



r(q- Ai)2(03(-q)03(q)). 



Because of the lattice symmetry C\ = C\ 



The scaling relation in Eq. 39 again remains the same 
provided we replace C with 3£a- As written, the integral 
for £a has to be done numerically over the momentum 
shell [shaded portion in Fig. 9]. 

Alternatively, an approximate analytical result of the 
integral can be derived by first doing a Taylor series ex- 
pansion with respect to in the denominators of Eqs. 57 
and then performing the integration. The result is shown 
in Appendix B. For small values of p, Eqs. 63, 66, 67 
can be expanded to obtain an approximate result for £a 
which is 



A-KsfpK 



1 



(59) 



(q • A2)2 + pi^(q • A3)2 + ^2(q . A2)2(q • A^f 



(57) 



This result shows that, as before, the shell integral di- 
verges as p 0. Once again, the scaling relation for p is 
identical to Eq. 39, with £ ^ 3£a- 

Finally, the evaluation of the constants A and xq 
which define the shell geometry in Fig. 9 are slightly in- 
volved, so we explicitly show how to evaluate them. We 
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take the real space basis vectors of the triangular lat- 
tice to be Ri ~ aox (along A3 in Fig. 8) and R2 = 
-(ao/2)x+ (V3ao/2)y (along Ai in Fig. 8). Then the 
basis vectors of the reciprocal lattice satifying the prop- 
erty Gi • Rj = 2TTSij are Gi = (27r/ao)x -|- (27r/-\/3ao)y 
and G2 = (47r/\/3ao)y. A is the distance from the Bril- 
louin zone center to the middle of one of the edges, and 
has length equal to half the length of the reciprocal ba- 
sis vector G2 (See Fig. 9). Thus A = 27r/(V3ao)- From 
the momentum shell diagram we find that the half- width 
of the Brillouin zone edge is xq = A/V^ = 27r/(3ao). 
Ai, A2, A3 are all equal to the lattice spacing ao which 
we have chosen to be of unit length. 



V. DISCUSSION: CRITICAL BEHAVIOR AND 
PHASE DIAGRAM 

Because the scaling relation for p is essentially identi- 
cal (up to numerical factors) for the calculations incorpo- 
rating the lattice symmetry and for the approximations 
used in Section IIIC, it is no surprise that the integra- 
tion of Eq. 39 gives results that are nearly identical to 
what we found previously. Fig. 10 illustrates typical re- 
sults for the square and triangular lattice case. As before, 
we find flows that either accumulate at the p = point, 
which we identify with the unbound dislocation phase, or 
flows that end at p > 0, which represent logarithmically 
bound dislocation pairs. As discussed above, this is a re- 
markable fixed point structure: there are no flows whose 
trajectory change discontinuously as the phase bound- 
ary is crossed; i.e., there are no relevant directions (in 
the KG sense) leading away from the p = flxed point. 
One can consider operators other than the ones we have 
discussed explicitly - cosines of higher order derivatives 
of the (/)'s, analogous to, for example, Eq. 41- but it is 
not hard to convince oneself that such operators will be 
more strongly irrelevant than those with which we have 
worked. Cosines of integrals of the ^'s are highly irrele- 
vant as well. 

In the absence of any such relevant directions, standard 
RG theory''' tells us that the free energy [F{£ 00)] is 
not singular, and one should not expect to find signals 
of a phase transition in thermodynamic quantities. De- 
spite this, there are some critical properties associated 
with the deconfinement transition as we have found it. 
For example, the screening length ^scr = loe'^ , defined 
by the scale £* at which a flow going to the unbound dis- 
location flxed point flrst strikes the p = axis (discussed 
in Section III C) , diverges as the confinement transition 
is approached. Eq. 39 suggests this divergence is a power 
law in (^Q — £,cr)t which may be confirmed by examining 
the flows in detail. The behavior of ^scr as a function 
of is illustrated in Fig. 11, where it is apparent that 
£scr diverges with a critical exponent —1/2 as 
from below. More explicitly, for flxed K this indicates 
^scr ~ \h — hc\~^^'^. The same power law is obtained for 
the triangular lattice. 



Another parameter of interest is the string tension 
which controls the effective logarithmic interaction be- 
tween vortices. Fig. 12 illustrates how this behaves as 
1^0 ^ Ccrl ^0. It is apparent that p vanishes quadrati- 
cally with this difference; or, equivalently, p ~ \h — hd'^. 
In principle the parameter p can be measured if one can 
create a vortex-antivortex pair at controlled locations in 
the system, and measure the force required to separate 
them. From the discussion in Section IIIB (see Eq. 37), 
it is clear that for a separation R, the force required will 
have the form (1 + C)/R, with C diverging when p = 0. 
The divergence indicates a change in behavior from a 
force that falls off with distance to one that remains flnite 
for any R, reflecting the presence of a string connecting 
the two vortices. We note that the diverging coeflicient 
C does not truly indicate a diverging force; even when 
p > 0, for small separations the force will be independent 
of separation, up to a characteristic crossover separation 
Rc ~ l/\/pK, where the cutoff in the denominator of 
Eq. 37 becomes apparent and the 1/R force sets in. Near 
the transition, the (l-|-C)/ii behavior only begins at very 
large R; the effective force between the vortices is thus 
never large. 

Finally, we discuss the reasoning leading to the pro- 
posed phase diagram illustrated in Fig. 1. For the pur- 
poses of this discussion, we consider K near its dual point 
{K w l/27r), which it should be noted for ft, = is well 
into the unbound vortex phase. For large h, the vortices 
will instead be linearly bound, even at this small value 
of K. As we know by now, with increasing l/h, there is 
a transition from an unbound to a logarithmically bound 
dislocation phase. This may alternatively be understood 
as a transition from linearly conflned to logarithmically 
bound vortices, as described in the last subsection. Since 
K is near the dual point {K « l/27r), the duality in 
the model tells us there must be another transition at 
large \/h as, Ec increases from zero, representing a tran- 
sition from unbound to logarithmically bound vortices. 
The dislocation deconflnement transition for which we 
developed the RG analysis is represented as the large Ec 
limit of the left transition line, and its dual is the large 
l/h limit of the right transition line. For vortices, we 
know that unbinding becomes increasingly difficult with 
increasing h, so that a smaller Ec is presumably required; 
this is why the phase boundaries move toward the axes 
as they approach the origin. We emphasize that the pre- 
cise behavior of these phase boundaries as l//i and E^ 
decrease is unclear, since our approach requires one (but 
not both) of these to be large in order to perform a con- 
trolled calculation. The diagram shown in Fig. 1 is the 
simplest that is consistent with what we have found in 
our calculations. As discussed above, the existence of 
the three phases is supported by numerical simulation 
studies^*. 



17 



VI. CONCLUSION 



VIII. APPENDIX B 



We have performed a renormalization group analysis 
to study vortex unbinding for the classical two dimen- 
sional XY model in a magnetic field on square and tri- 
angular lattices. Our analysis shows that if one starts 
at high temperature and large magnetic field, vortices in 
the model unbind as the field is lowered in a two-step 
process: first strings of overturned spins proliferate, and 
then vortices unbind. The proliferated string phase may 
be understood as one in which the vortices arc logarith- 
mically bound, whereas in the high field phase they are 
linearly confined. This transition has an alternate de- 
scription in terms of domain walls and screw dislocations, 
at which the dislocations deconfine. The vortex decon- 
finement transition was shown to be dual to this. Both 
transitions arc remarkably continuous, but they arc not 
of the Kosterlitz-Thouless type. Finally, the unbound 
vortex fixed point was shown to contain a set of nodal 
lines which are lines of zero energy modes, reflecting the 
symmetry of the lattice. 



VII. APPENDIX A 

In this Appendix we display the explicit form the shell 
integral in Eq. 34. Defining the parameters A = 1 + B, 
B = Kp + A^^, one finds 



2KA? 



1 



X 



i?+ - 1 1 

— arctan — 

R+ /R^ 

i?_ - 1 1 

arctan — ;= 

R- ./Ri 



+ 



B 



1 B 

arctan \ — 

/AB V A 



(60) 



In Eq. 60, the parameters R± are given by 

Eq. 60 can be expanded for small p to yield the result 



2iCA2 I TT 



+ ■ 



1 



1 + ^2A2 
1 



Ae(l + C'A2) 



arctan 



AC 



(1 + A2^2)3/2 



(61) 



For small values of p, the first term in Eq. 61 dominates 
the RG flows. 



In this Appendix we show the approximate analytical 
results for the following shell integral (See Eqs. 58) with 
p, = pK 



C 72 

Jsholl l^"" 



(q- As 



(q • Ai)2 + ^(q . A2)2 +ei^■ Ai)2(q • As)^ 



(62) 



To evaluate this integral we expand the denominator for 
small ^2. The integral now looks as. 



£ A {l^, Odi = £X il^)d£ + -T^A (m)'^^, (63) 



t2 d ^(b). 



where, 



jC''2^de = K 



L^^dl = K 



f ^'q 

ysheii(27r)2 (q 



dp, 



(q- A2)2 



Ai)2+^(q. A2)2' 
d\ (q-Ai)2(q. A2)2 



/shell (27r)2 (q.Ai)2+Mq-A2)2- 
We now display the results of the above integrals. 



l + 4p 



4XV3A2 / 3A2 + 2V3Axo + xg(l + 4p) 
~ (1 + 4/^)2 \^3A2-2\/3Axo+a;^(l+4/z)^ 

2KV3A'^(l-4p) ( aVsAxoJp 
H , ■ — arctan ' ^ 

2KAxo 2V3KA'^ 
+ , , + 



3A2-a;2(l+4M) 



X In 



1 + P ' {l + P^ 

^3A2(1 + p) + 2V3Axoip - 1) + a;2(l + //)' 

3A2(1 + p)^ 2VSAxo{p - 1) + a;2(l + p) 



2KV^A^(1-p) ( 4V3Aa;oVM 

H ^rrz — — TTT^ arctan ' 



\{l + p){3A^-xl)^ 

2KAxq 
4-H^ 

4^/3A2 2,A^p + 2^/2,Axop + xl{4: + p)^ 



2-K 



dp 



+ 



8j^AA2xg _ 24JfA3A2(16/z2 - 2,2p + 2,)xo 

~3(4/i+ 1)2 {Ap + ly 

24V3Js:A'^A2(32^2 _ 2Qp + 1) 
(4/x + 1)5 



(64) 
(65) 



{A + pY \^ 3A2y^ - 2 VSAxoAi + xl{A + p) ) 
2KV3A\p-A) ( AVSAx,^ \ 

vni^+pr ^'''''''[sA^p-xiiA+p) ■ 
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X In 



+ 




+ 



3A2 + 2V3Aa;o + a;§(4/x 
^3A2 - 2\/3Aa;o + a;2(4/i + 1) 

{4n + 1)4[[3A2 + xliAfx + - 12A^xl] 
( 2,^/2,K^^^?{l'j2^l^ - 560/i2 + 100/i - 1)' 



X arctan 



3A2-a;2(4M+l) 



36Ji:A5A2[3A2 + xljAn - l)][16/x2 - 24/x + 1]xq 
(4m + 1)4 [9A4 + 6A2a;2(4M - 1) + a;^(4;U + 1)2] 

6(/i + l)2 

12V3KA4A2(2y^2 



5/z + l) 



In 



+ 



^ 3A2(l + /^) + 2^/3A(/z- l)a;o + a;g(M+l) \ 
3A2(1 + /i) - 2\/3A(/x - l)xo + xlin +1)) 

288if A5A2(3A2 + xl)xon{^i - 1) 
(m + 1)4[[3A2(1 + + (m + l):r:g]2 - 12A2x2(m - 1)2] 

^12V3i4:A4A2(3/x3 - 35/x2 + 25At- 1)^ 



X arctan 



4^/3A^a;o 



(l + A*)(3A2-ar2)y 
{l44iirA^A2(3A2 - a;g)a;o(Ai - 1)(m^ - 6/i + 1)} 

{(/^ 



1) 



9A4 - 6A2 



x;^(/x + 1)2 + 18AV 



36A2a;§/i - QK^xln^ + 9A'^/i2 

144ifA3A2(/x-4)xo 



2ifAA2xg 



+ 



3(m + 4)2 
48V3ii'A4A2(/i2 



(m + 4)^ 
10/i + 8) 



(a^ + 4)5 
3A2/i + 2\/3A/ia;o + .x^(m + 4) ' 
I 3A2m - 2 VSA/zxo + + 4) 



+ 



+ 



1152ifA5A2/z(At - A)xl 



+ 4)4[[3A2m +{^Ji + 4)a;2]2 - 12A2/z2a;2] 
^ 3\/3JrA4A2(3/x^ - 140/^2 + 400^ - 64) ^ 



X arctan 



VJ^(P + 4)5 
/ 4V3A^.T o \ 



\^3A2M-a;2(/x + 4)^ 
{36Js:A5 A2 [3A2m + 2;2 (4 - /i)].To - 24/x + 16] } 

+ 4)4 9AV^ + 2AK^^ixl - <q1^p?xI 
4(m + 4)2]}"\ 



(67) 
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FIG. 11. 1/^scr VS. ~ icr\ for squaTc lattice with full 

symmetry where Iscr is the screening length. Iscr diverges 
with a critical exponent —1/2 as ~^ Ccr from below. Pa- 
rameters are same as in Fig. 10. 
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FIG. 10. RG flows for the scaling relation in Eq. 39 for 
(a) square lattice and (b) triangular lattice with full lattice 
symmetry. K = 0.01 and initial po = 1/40. Left vertical axis 
is a fixed line, and the thick line is a separatrix between flows 

that reach p = for finite £ and flow to the origin as ^ cx), 
and those that have p > at the end of their flow. Thick lines 
are separatrices with critical values of (a) = 7.755 and (b) 
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FIG. 12. Normalized string tension p/po as a function of 
(^0 ~ Ccr)^ for the square lattice with full lattice symmetry. 



p vanishes quadratically as 
are same as in Fig. 10. 
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